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Abstract. Solvable vertex models in a ferromagnetic regime give rise to soli- 
ton cellular automata at q = 0. By means of the crystal base theory, we study 
a class of such automata associated with the quantum afBne algebra C/q(gn) 



for non exceptional series Qn = A 



(2) 



i(2) 



dI^' and D. 



(2) 



They 



possess a commuting family of time evolutions and solitons labeled by crystals 
of the smaller algebra Uq{g„—i). Two-soliton scattering rule is identified with 
the combinatorial R of l/q(g„-i)-crystals, and the multi-soliton scattering is 
shown to factorize into the two-body ones. 



1. Introduction 



1.1. Background. In early 90's the crystal base theory explained |KMN| the 
appearance of affine Lie algebra characters | DJKMO{ in Baxter's corner transfer 
matrix method ||^]. The success concerned solvable vertex models in an antiferro- 
magnetic regime. Typically for the six- vertex model corresponding to Uq{A^^"'), the 
spins in this regime are alternating at q — 0. On the other hand the ferromagnetic 
regime, where the distant spins are all equal, had been left intact until rather re- 
cently. A renaissance arose with the discovery |HKTl| , |FOY| , |HHIKTT| , [HKT2| 
that the ferromagnetic vertex models at q = turn out to be soliton cellular au- 
tomata pTSl , As a result the ultradiscrete soliton equation |rTMS] acquired 
an origin also as a quantum integrable system at q = 0. 



In | HKT1[ a class of soliton cellular automata were constructed associated 



A 



(2) 



(2) 



^(1) 



The solitons therein possess internal degrees of freedom labeled by 



with [/^(g„)-crystals for non exceptional series g 
and D^^l^ 

crystals of the smaller algebra U^(Qn-i)- 

The purpose of this paper is to study the scattering of the solitons in these 
automata. The two-body scattering is characterized by the data: 

• Exchange of internal degrees of freedom. 

• Phase shift. 



In |jFOY| it was shown for g„ = An '' that they fit completely with the combinatorial 
R and the associated energy function for Uq{A^j^}_^)-CTystals. Moreover, the multi- 
soliton scattering factorizes into two-body ones. See also [rNS|. 
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The main result of this paper is Theorem 2.7, which extends such results to 



the 2n automaton |HKT1| for all the non exceptional g„ mentioned above. In these 



automata local states take values in the crystal of the vector representation. We will 
also consider the generalized automata corresponding to the analogue of symmetric 
tensor representations of U^{gn)- Except for the phase part, such a generalization 



has been studied for An^ case in |HHIKTT[ . Our proof of Theorem |2.7| contains 



a new ingredient, an operator T\f, which transforms a soliton into another. To 
illustrate the idea, let us prove the An ^ case in the next subsection. 

1.2. case. Let us quickly review necessary notations. For U^lAn"^)- 

crystal Bi, see Appendix First z; is a map sending an element (xi, a;2, • • • , a:„) of 

the C/^(A|^^li)-crystal Bi to |n+l [ ^"""g)- • •^[i]®''' (g>\2f''' in the ?7^(Ai^^)-crystal 



Bf . It is known [ KN that z; is an isomorphism of ?7g(A„_i)-crystals in the sense 
that ii ■ e-i — e^+i • ii and z; • fi = /i+i • z; for z = 1, 2, • • • , n — 1. A one-soliton state 
of our cellular automaton is described as 

••■«)|T|(8)---®|T|®z/(&)(8)|T|«'---(8)|T|(8)... 

for some b £ Bi. Such a soliton is said to have length I. Under the time evolution 
Tr (see Section 2T), the length I soliton moves rightward by min(Z,r) steps. One 
can now define the phase 7 by counting the position ^(t) at ti me t from the right 
end as — min(Z, r)t + 7 imder the time evolution Tr. See Section for details. To 
the one-soliton z/(6) {b e Bi) of phase 7, we associate an element z'^b of Aff{Bi) := 
{z'^blbe Bi,dGZ}. 

Consider a two-soliton state z'^^bi (g) € Aif(i3/J (g) Aff(i3;J {h > ^2,71 > 

72) illustrated as 

....ti^{bi)....ii^{b2) 

If r is greater than I2, moves faster than ii2{b2) under Tr. After sufficiently 

long time, we have 

HAb'2)-MAb'i).:. 

which contains solitons of the same length U (z — 1, 2) as the initial state. Let 7,' 
be the phase of the soliton iii{b[). Now the theorem is 

Theorem. The scattering rule: 

z^^bi ® z'>^b2 I — > z'^^b'^ (g) z'^'^b'^ 

is given by the combinatorial R : 

R : z^i6i ® z^^b2 I — z'^-'+^b'^ ® z'^i-'^fo'i, 

where 5 — H [bi ® b2) is the energy function. (Therefore 7^ = 71 — <5, 72 = "f2 + S.) 



The proof goes as follows. First note that Tr and R both commute with the 
actions of e^, for z = 1, • • • , n — 1 on Bi^ ®Bi^ or Bi^ ®Bi^ (t/q(A„_i)-invariance). 
Thus we are allowed to reduce the proof to the set of highest weight elements, i.e., 
elements killed by (z = 1, • • • , n — 1). The set of such elements in Bi^ ® Bi^ is 
denoted by H{Bi^ ® Bi^). If li > I2, we have 

HiBi^ (E)Bi^) = {ui^ (g) (a;i, 2:2, 0, • • • , 0) | + a;2 = ^2}, 

H{Bi2 (gSii) = (g (yi, 2/2,0, • ■ • ,0) 1 yi +y2 = ^i,y2 < h}, 
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where ui stands for 0, • ■ • , 0) £ i?;. It is known that 
R : Aff{Bi,)(g, AfF{Bi,) AS{Bi,) d) AS{Bi,), 
z'f'ui, (dz^^xi,X2,0,--- ,0) I — > z'<^+"ui,(g)z'"-"{xi+h-l2,X2,0r-- ,0), 



where H — X1 + I2 under suitable normahzation of the energy function. (See (A.l).) 



For the proof we make use of another operator Ti, constructed from the {An )- 
crystal , which corresponds to the degree 2 antisymmetric tensor representation. 
Tt| is to put the highest element (2) G B[^ on the left and to carry it throu gh to the 



right via the isomorphism (g) Si (g) ■ • ■ (g) Bi ~ i3i g) • ■ • (g i3i (g . (See | NY for 
the algorithm to compute g) Bi ~ Bi g) Bt|.) It acts on the space of states and 
has the following properties: 

• commutes with T^. 

• On a one-soliton state of length I, acts as 

rri 11 n ^^^ / ^''l^, 0, 0, • • • , 0)) if 2:2 = 

If one cuts the figure given below in the middle, the left and the right pieces become 
examples of 0:2 = and 2:2 = 1 cases in the above, respectively. 

Now suppose our theorem is valid for a particular highest weight element 
z'^^ui^ (g z'^'^{xi,X2,0, • ■ • ,0) (x2 > 0), i.e., 

Tl{z'^^ui,®z'-<^{xi,X2,0r-- ,0)) 

= z^^+^i+'^uj, ® + ^1 - Z2, a;2, 0, • ■ • , 0) 

for sufficiently large t. Apply on both sides and use the commutativity T^T^ — 
TrT^. We get 

Tl{z^^ui, (g 0-^^-1 (m + 1, a;2 - 1, 0, • • • , 0)) 

= z-^^+^i+'^w,., +Zi -/2 + 1,2:2 - 1,0,--- ,0). 

This shows our theorem is also valid for z^^ui^ g) z''^^^{xi + 1, 2:2 — 1, 0, ■ • ■ ,0). 
The following figure explains the action of Ti, on a two-soliton state (It transfers 
the upper line into the lower line): 

T^(z'^iw3 02-^^(1,1,0,- •■ ,0)) = 2^1^3 0)2-^=^-1(2, 0,0,- •■ ,0). 



2221113211 
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22222223333 
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One should notice that the above procedure is invertible. We thus reduced 
the theorem to the simplest case z''^ui^ (g z^'^ui^, where the scattering rule can be 
verified directly. 

1.3. Plan of the paper. In the rest of the paper we simply write g to mean 
g„, and let g(= g^) = Cn,Cn, Bn,Cn, Dn and Bn denote the classical parts of 
— A'^J^-n ^^211, J ■^ri\ C'n\ Dn^ and f^+i, respectively. As for the crystals we 
will use the same symbols Bi and Aff{Bi) either for q or g„_i as the distinction 
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is clear from the context. Bi and AfF{Bi) will be called ?7^(g) and C/q(0)-crystals, 
respectively. 



Our proof for the listed g is basically similar to Section 1.2. We start with 
dividing the cases into three types: 

Type I 

Type II dn^ 
Type III 4l\^i+i 

In each Type the decomposition of Bi or B\^ into J7g(g)-crystals and the parameter- 
ization of the highest weight elements of Bi-^ ® Bi^ are the same. 

Bi = B{lAi) for Type I, 

= B{lAi) © B{il - 2)Ai) ® • ■ • © (B(Ai) or B{0)) for Type IL 

= B{lAi) © B{{1 - l)Ai) © • • • © 5(0) for Type III. 

= B(A2) © B{Q) for Type I, 

= B{A2) for Type II, 

= B(A2)©B(Ai)©B(0) for Type III. 

For Type I and II the proof proceeds just as in the An'' case. The only difference 
is that Tt| is not invertible. But it can be remedied by virtue of the fact that the 
weights in the inverse image are distinct. For Type III the proof is reduced to Type 
II by using some embedding of [/^ (g)-crystals into J7q(ci^^)-crystals. 

In Section || we recall the basic feature of the g automaton. Soliton states and 
their label in terms of tensor products of AfF{Bi) are explained. Our main result 
is Theorem In Section || we give a proof along the idea explained in Section 
tioiTl 



1.2. In Section M typical examples of the scattering of solitons are presented. In 
Section ^ we consider the generalization where the automaton states are taken from 
• ■ ■(S)Bmj 'X'Bmj+i ® • ■ ■ • We show that the inhomogeneity of {rrij} affects the soliton 
scattering only via a certain phase shift. Appendix H is a summary of basic facts 
on crystals and combinatorial R . Appendix ^ and ^contain some information on 
the crystals Bi and B\f, respectively. 

We close the section with a remark. Solitons in our automata with an equal 
velocity can form a bound state (breather) in general. It propagates stably if 
isolated, and exhibits curious reactions with other solitons or breathers. In this 
paper we shall only study the scattering of solitons with distinct velocities, and 
leave the study of bound states as a future problem. 



2. Soliton cellular automata 

- 4(2) 



2.1. States and time evolutions. For g = A^^^_i(n > 3), Al,^„'(n > 2), B^^^(n > 

3),D'h \n > 4),D';^l^{n > 2), let Bi be the C/^(g)-crystal in [[KKMl] . For g 

Ci^\n > 2) our Bi is the one in ||HKK0T] , |HK0T1| . We have hsted the param- 
eterization of Bi as a set in Appendix y. For any g we write 

(2.1) ui = {l,0,... ,0)eBi 

in the notation therein. 
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Consider the crystal Bf^ for sufRciently large N. The elements of Bf^ we 
have in mind are of the following form: 

■ ■ ■ (S)[l]<S) ■ ■ ■ <E)[l\(E) bi (g) ■ ■ ■ (g) bk <^[l\(E) ■ ■ ■ ■ ■ ■ , 

where bi, . . . ,bk £ -Bi and [T] — ui. Namely, relatively few elements are non [T] , 
and almost all are [T]. See Appendix ^ for the definition of the symbol [a]. In the 
assertions below, we embed, if necessary, Bf^ into Bf^ (N < N') by 

ci (g) ■ • ■ Cat i-^ ci ■ • ■ (g) Cat (g) [T] (8) ■ • ■ ® [T] . 

^ ^ 

N'-N 

Let _B; (g) -Bi — i?i Bi be an isomorphism explained in Appendix ^ 

Lemma 2.1. By iterating Bi (g Bi ^ Bi (g Bi, we consider a map 

Bi®Bi®---®Bi^Bi®---®Bi®Bi 
ui®bi® ■ ■ ■ ®bN 1-^ bi® ■ ■ ■ ®bN g)b, 

where we assume bj = |T] for N' < j < N . Then there exists an integer Nq such 
that b = ui for N - N' > Nq. 



Taking sufficiently large N — N' such that the above lemma holds, we define a 

N- 



map Ti : Bf^ — > Bf" hy bi (g ■ ■ ■ g) bN ^ bi (» ■ ■ ■ g) b 



Lemma 2.2. For a fixed element of Bf^ as in Lemma 2.1, there exists an 
integer Iq such that Ti — Ti^ for any I > Iq. 



Both lemmas can be verified directly by using the result in | HKOTl , HKOT2 1 . 
In particular m = 1 case of the fact 

(2.2) ui (g) Um — Um g) ui Under Bi eg) Bm — Bm ® Bi for any /, m, 
is relevant here. 

An element of Bf^ having the property described in the beginning of this 



subsection will be called a state. Lemma 2.1 and 2.2 enable us to define an operator 



T = lim;^oo Ti on the space of states. Application of T induces a transition of state. 
Thus it can be regarded as a certain dynamical system, in which T plays the role 
of 'time evolution'. We call it g automaton. By the same reason, Ti may also be 
viewed as another time evolution. (In this paper, time evolution means the one by 
T unless otherwise stated.) 



Under the condition stated after Lemma 2.1, the sequence 6i g) • ■ • g) is 
determined uniquely from 6i g) ■ • • gi^Af and b = ui, since the isomorphism Bi®Bi ~ 
Bi g) Bi is bijective. Hence, the time evolutions T; {I > 1) are invertible. 

2.2. Conservation lavirs. Fix sufficiently large N and consider a composition 
of the combinatorial i?'s 

Ui ^Rnn+1--- R23R12 : AfFiBi) ® Aff(Bi)®^ AfF{Bif^ ® Aff(Bz). 
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Here Ra+i signifies that the R acts on the i-th and (i + l)-th components of the 
tensor product. Applying TZi to an element ui®p {p = bi ® ■ ■ ■ ® 6jv), we have 



N 



where 6*^°^ — ui and b^^^ (1 < j < N) is defined by 



Bi®Bi 



Bi ~ Bi 



' Bx ®Bi 



] 3 

ui®hi® ■■■ ®bj ^ hi® ■■■ ®b4 ® b^^^ . 



Z°bN 



z"Nbj, 



Since H{ui ® ui) — 2i; ^ (see (AJ)), the quantity Ei{p) grows linearly with 
N. To avoid it we introduce a regularized version: 



AT 



Blip) ^ - J2{Hj - H{ui ® Ul)) = Ei{p) + NH{ui ® ui), 



which becomes independent of N for iV sufficiently large. 

The Yang-Baxter equation in Proposition [A.1| implies that as maps from 
AS{Bi) ® AS{Bi,) ® Aff(Bi)®^ to AS{Bi)®^ ® Aff{Bi,) ® AfF{Bi), the equality 



(2.3) 



TZi'TZiRu = Rn+1 N+2T^iT^i' 



is valid. Here TZi (resp. 7^^0 acts on AfF{Bi) ® AfF{Bi)'»'^ (resp. Aff(B/-) 
Aff(_Bi)®^), and acts on the other component as identity. 



Proposition 2.3. For an element p e Bf^, we have 

(1) TiTi'{p)^TvTi{p). 

(2) Ei{Tv{p)) = Ei{p). In particular, Ei{T{p)) = Ei{p). 



Proof. We follow the argument in Section 3.2 of |FOY|. Consider the element 
z°ui®z°ui' ®p G Aff{Bi) ® Aff{Bi,) ® Aff{Bi)®^ . Apply the both sides of to 
it. Graphically one gets {6 = H{ui ® ui')) 
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Tiip) 



Tvip) 




TiT,,{p) 

Since Ei{p) differs Ei{p) just by a constant, the assertions (1) and (2) follow directly 
from this diagram. I 

In addition to Ei, there is another conserved quantity resulting from the 
invariance explained in the beginning of Sec tion ^. However we do not detail it in 
this paper. For An '' case, see Section 4.2 of |FOY|. 



2.3. Solitons. In order to describe solitons in the g automaton, we introduce 
an injection 

ti : Bi for {/;(0„_i) ^ (Bi for {/^(g))®' 
for each / G Z>i as follows. 

3 = ^2n-i ■ Foi' ^ = i^i^ ■ ■ ■,Xn-i,Xn~i, . . . in L/^ (0„_ i )-crystal Bi, 



rzz]<»x„^i I iQ 

) | n| ® 



g = A^,; : For 6 = {xi, . . . ,a;„_i,x„„i, . . . in [/^(g„_i)-crystal Bi, we define 
s{b) = ^'(fe) = [{1 - s(^))/2]. 

If I — s{b) is odd. 



»i(6) = 0«)|T] _ 
otherwise 



(gixi I |(g)s'(f)) 



I iQ9a;„--i I ieS2;„_i I i(»xi I 1 

\n\ (g)[n\ • ■ • |_2j (E>[l] 



■10X1 \~~\^s'{b) 



= Si^' : For 6 = (xi, . . . , a;„_i, xo, ^„-i, ■ ■ ■ ,xi) in [/^(g„_i)-crystal B;, 

(gixo I \0x„-i 



n{b)=2 



nnQ9x„_i „r--\'»xo I |Q 



)0 



g = C,l^^ : For 6 = (xi, . . . , .Xn-i, a::„_i, . . . in [/^(g„_i)-crystal B;, we define 
g = Z?!^' : For 6 = {xi, . . . , a;„_i, x„_i, . . . in C/^(g„_i)-crystal Bi, 

,,, r^®Xi I \0X„-^i I |l8x„_i I llglXi 
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D 



(2) 
n+1 



: For b 



ci) in [/^(g„_i)-crystal Bi, we 



define s{h) = a;o + ELi (a^* +3^0, s'(^) = [(^ - s(fe))/2]. 

U] ®---®[^ (»[0] 



If / — s(6) is odd, 
otherwise 

n«is'(b) 



)0 ®g 



The ii here was originally introduced in |HKTl| , where the following was 
shown. 

Lemma 2.4. For i e {1, 2, . . . , n — 1}, the following diagram is commutative: 



Ug{gn-i)-crystal Bi 

-'A 

U'{Sn-i)-crystal U {0} 



U'{q)- crystal Bf 



ei+i 



U:j{Q)-crystal U {0}, 



where we set ii{0) = 0. T/ie same relation holds between ft and fi+i- 

The Kashiwara operators of Ug{Qn-i) and J7q(0)-crystals should not be con- 
fused although we use the same notation. The proof is due to the explicit rules for Ci 



KKM, HKKOTI and the embedding of ?7'(fl„_i)-crystal Bi into f/'(0„_i)-crystal Bf 



as f7g(fl„_;^)-crystals |KN |. 
A state of the form 

(2.4) 



is called an m-soliton state of length li,l2, ■ 
figuration such as 

••■(g)l(8)l(g)z;(5)(8)l®l(g 



Jrn- Here .. 



denotes a local con- 



for some b E Bi 



surrounded by sufficiently many I's. In view of (2.6) below, we assume that there 
are at least k I's (i.e., 1®'*) between [k] and [h+i]- (Here and in what follows, we 
often abbreviate [a] S i?i to a.) In the sequel we set 



(2.5) 



1 otherwise. 



By combining the results in pKOTlj , |HKOT2| , [HKTl[ it is not difficult to show 

Lemma 2.5. Let p be a one-soliton state of length I, then 

(1) The k-th conserved quantity of p is given by Ek{p) = <;min(k,l). 

(2) The state Tk{p) is obtained by the rightward shift by mhi(k,l) lattice steps. 

For example (2) is shown by establishing 
(2.6) uu ® nib) ® l^^C^'O ~ i®™n(fe^O n(b) ® Uk, 



for any b e [/^(0„_i)-crystal Bi under the isomorphism Bk (8) B\ 

Bf 



®/+niin(fc,/) 



q 
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For any state p, define the numbers Ni ^ Ni{p) = 1, 2, . . . ) by 



k>l 

Ni = {-Ei^i+2Ei-Ei+i)/c;. 



By Lemma 2.5, we have 



Proposition 2.6. For the m-soliton state ({2.4 ), Ni is the number of solitons 
of length I, i.e., Ni = ft{j | Ij = I}. 

This proposition implies the stabihty of solitons, since the numbers Ei{p), hence 
Ni{p), are conserved. 

2.4. Scattering of solitons. We introduce a labeling of solitons of length I 
using Aff{Bi) for the lower rank algebra [/^(0„_i). Suppose there is a soliton of 
length I at time t 

(2.7) • • • (g) 1 ® ii{b) ®l(g)--- = -- -(g)l(8)j/i(g)z/2«'---®i'/®l(8)---, 

where vi,... ,1^1 G Bi and b is an element of [/^(g„_i)-crystal Bi. Say it is at 
position 7(t), if I'l is in the 7(t)-th tensor component of Bf^ counting from the 
right end. From Lemma ^.5| (2), the position j{t) under the time evolution Tk 
is given by 7(t) = — min(fc, l)t + 7 unless it interacts with other solitons. Here 



min(/c, I) is the velocity and 7 is the phase. To such a soliton (2.7) we associate the 
following element in Aff{Bi). 

(2.8) Hq^A^^Id^^I,, 

(2.9) z^n if0 = 4^j,i5i^|i, 

(2.10) z^-'-'b if0 = 4^),i?i'|i, 1.1 = 0, 

where the symbol (f> € Bi has been introduced in Appendix ^ 
Now consider an m-soliton state illustrated as below. 



As in (2.4) we assume that solitons are enough separated so that there are at least 
li I's between [li] and [h+i]. To such a state we associate 

(2.11) zTifci ® z'^^b2 • • • (g z'f-^b^ e Aff{Bi,) (g) AfF{Bi,) ® • • • ® Aff{BiJ, 

where each element z'^^bj G Aff{Bi^) is specified as in the previous paragraph. 
Suppose h > I2 > ■ ■ ■ > Im- Since longer solitons move faster, we can expect that 
the state turns out to be 



after sufhciently many time evolutions. (This claim is due to Corollary 2.S 
below.) We represent such a scattering process as 

(2.12) ^ z''''-b'^(g)---(g)z'''^b'2(g>z'''^b[. 

The main result of this paper is the following. 
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Theorem 2.7. The two-body scattering of solitons ii{bi) and ik{b2) (I > k) 
under the time evolution Tj. (r > k) is described by the combinatorial R : 

Aff{Bi) ® AS(Bk) - AfF{Bk) ® Aff{Bi) 



^61 



with S = H{bi 62). 



The statement on the phase shift is vahd only for the energy function normaUzed 
as (|A.l|) . Such values wiU be given in (2) of PROPOSITIONS |t], js^ and Theorem 



3.10 



From Theorem 2.7 one can easily establish the followin g with the help of 
Proposition 2.3 and Lemma 2.5 by the same argument as in [ FOY | . 



Corollary 2.8. Scattering of solitons is factorized into two-body scattering. 

the iso 



Namely, (2.1i) is valid under the isomorphism of crystals 
Aff{Bi,)®-- 



AS{Bi_ 



AS{Bi,). 



The map decomposes into the combinatorial i?'s corresponding to pairwise 
transpositions of the components. It is independent of the order of the transposi- 
tions due to the Yang-Baxter equation (Proposition A.l ). 



Although we have assumed Zi > ■ • ■ > in Corollary 2.8, some part of the 



results can be generalized to the situation where there exist solitons with the same 
length. See the second remark after Theorem 4.6 in |FOY|. However a thorough 
study of such cases will necessarily involve the bound states as mentioned in the 
end of Section 1.3. 



3. Proof of Theorem 

3.1. Preliminaries. First we recall an important consequence derived from 
representation theory. Note that U'^{q) contains C/q(0„_i) as a subalgebra. This 
fact can be translated into the language of crystals, and combined with Lemma 
2.4, it leads to the property that on an m-soliton state 



Tk commutes with e^, (i = 1, 2, • • ■ , n - 1) on Aff(i?; J Aff{Bi^). 

Here we are identifying the m-soliton state with an element from Aff{Bi-^) ® • • • (Xi 
Aff{Bi^ ) as explained in Section O. The action of Tk on Aif(i?; J ®---® AS{Bi^ ) 
is defined through this identification. The actions of Ci/fi (i = 1,2,--- ,n — 1) 
are calculated according to the rule explained in Appendix [A|. By the actions, 
the power of z in an element of AfF{Bi . ) is unaffected. We c all this property 
[/q(5„_]^)-invariance. This reduces the proof of Theorem 2.7 to the subset of 
Aff{Bi) (g) Aff{Bk) corresponding to the highest weight elements with respect to 

Uqidn-l)- 

So we give the data on the highest weight elements of Bi (g) Bk and their image 
under the combinatorial R. For elements in Bi of a particular form, we employ the 
following notation: 

(3.1) (xi,a;2,xi) = (a:i,a;2,0, ••• ,0,xi), (a;i) = (a;i, 0, • • • , 0). 

For Type I case, we have 

Proposition 3.1 ( pKOTl] , |HKOT2[ ). Assume I > k. 
(1) The set of highest weight elements of Bi ® Bk is given by {{I) ® {d,c^b) \ 
b + c + d^ k}. 
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(2) The combinatorial R acts on the set of highest weight elements as 
(3.2) R : (g) c, b) ^ 0"'+^(fc) (g) z'^-^il -b^c,c, b), 

where H = 2k — 2b — c. 
For Type II case, we have 



Proposition 3.2 ( pKOTlj ). Assume I > k. 

(1) The set of highest weight elements of Bi ® -Bfc is given by {(/) ® {d,c,b) \ 
f <IJ = 1 {mod 2),b + c + d<k,b + c + d = k {mod 2)J>b + c}. 

For a highest weight element (/) (g) {d,c,b), define nonnegative integers a,e by I = 
2e + f,k^2a + b + c + d. 

(2) The combinatorial R acts on the set of highest weight elements as follows: 

(i) Ifa>e, 

R : z"\f) ® z"''{d, c, b) z"'+^(fc - 2e) «) z"'-"{d + I - k - y,c,b ~ y), 

where y — min(^ — k,{b — H — k + a~e + {d — 

(ii) If a < e and d — b<e — a<l~k, 

R:z"'{f)(E)z"''{d, c, b)^ 
z"''+"{k - 2a) (g> z"'~"{d + l-k-e + a- z,c,b + e-a- z), 

where z — min(fe — d + e — a,l — k — e + a), H = k. 

(iii) Otherwise, 

R : z"\f)(g)z"''{d,c,b) ^ 

z"''+"{k - 2e + 2w) ® z"'-"{d + l-k-w,c,b + w), 

where w = m\n{l — fc, (2e — 2a — d + 6)+), H — k + max(e — a — I + 
k,d — b — e + a). 

Here and in what foUows we use the notation 

{x)± — max(0, ±.t). 



Lemma 3.3 ( pKOTll , |HKQT2[ ). Under the isomorphism Bi®Bi ~ Bi®Bi 
we have 

^ \0®a,o,o) j = o, 

^ 10® (0,^,0) j = l. 

This is vahd for all Types I, II and III, and it is essentially the same with the 
An'' case. 

We also give all the images of the isomorphism B^ ® Bi — > Bi ® B\^ that will 
be needed in the next subsection. Let the symbol 
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signify (^) ® [7] [y] ® (^/) under the isomorphism, li i ^ jihj = 1, 2, ■ • • ,n), 
we easify obtain 



from the weight reason. We hst up the other cases below. 



Common for Type I and II: 



3 3 
For Type I onfy: 



1 1 
For Type II only: 



3.2. Type I and II. Let p G Bf^ be a state. For an element v G B\^, define 
a map : Sf ^ Sf ^ by 



(3.3) 



® B^ 



where v' is uniquely determined from v and p although it is not necessary for 
defining [v] . 

For Type I case, we have 
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Lemma 3.4. Consider the one-soliton state identified with z'^{d,c,b) as an el- 
ement of AS{Bk). Then we have 



z"'-'^(d+l,c.h-l) ifh>0, 
{z"'{d,c,b)) = { z"-i(d+l,c-l,0) ifb = 0,c>0, 
z"'{d,0,0) ifb = c = 0. 



For Type II case, we have 

Lemma 3.5. Consider the one-soliton state identified with z'^{d,c,b) as an el- 
ement of AfF{Bk)- Define a by k = 2a -\- b -\- c -\- d. Then we have 
lfa = b = c = 0, then 



{z^{d,c,b))=z"'{k,0,0), 



else 
2 



l,c,5- 1) 



if bd > 0, 



{z"'{d,c,b)) = i (6- l)_,c- 1,(6- 1)+) ifbd=0,c>0, 



^(d+(6-2)_,0,(fe-2)+) 



ifbd = c = 0. 



{z"^{d,c,b)) 



z"'-\d+l,c,b+l) ifa>0, 

l,c,6+ 1) ifa = 0,d>0, 

z™(0,c- 1,6+ 1) ifa = d = 0,c>0, 

z'"(0,0,6) ifa = d = c = 0. 



Instead of proving these lemmas, we give examples below. 
Example 3.6. (i) Type I (A; = 4) 



Thus we have 



(0'"(l,l,2))=z'"-^(2,l,l) 



(ii) Type II (fc = 6,a=l) 



1 1 



Thus we have 



(^'"(0,l,3)) = ^™-i(0,0,2). 
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(iii) Type II (fc = 6, a = 2) 



Thus we have 



(z™(2,0,0)) = z'"-i(4,0,0). 



In what foUows, we understand T[, = T], [(2)]. Let 
Pi = «fci(&i) 

P2 = «fci(&l) «fc2(&2) 



be the one and two-sohton states which are identified with z"^^bi G Aff{Bki) and 
&i(8)z™262 e Air(SfcJ(g)Aff(BfcJ, respectively. Suppose that (2) ^Pi ^ T^{pi 



v' as in (3.3). Then Tif(pi) is another one-sohton state identifi ed w ith some element 
rt,(z™i&i) e AfF{Bki), which is exphcitly given in Lemmas and 3^ for 61 of 



the form Similarly, T[^{p2) is a two-soliton state identified with T[^{z"^'-bi) (g) 



Tt,[i;'](z™2 62) e Aff(BfeJ Aff(Bfe.J. When 61 (g) 62 is a highest weight element, we 
will only encounter v' = (2) or (^) . In fact one always has v' — (2) for Type I (and 
An ^). On the other hand for Type II, v' — (^) happens if and only if 61 — (/) with 
/ < fci - 2, and v' = (2) otherwise. See Proposition (1) for the form of the 
highest weight elements and Example 3.6 (iii) for a situation v' — (^). Thus the 



results on Tti[(2)] -^J©] Lemmas 3^ and 3^ are enough to calculate the 



action of iVfi)] ^"^V highest two-soliton state z"^^bi (8) z"^^b2- Later in the proof 

[either for {ki, = {I, k) 



of Lemma |3.8| , we will extensively use Lemmas and | 
or (fc, I). 

Let us prepare a few more facts. 



Proposition 3.7. 



Proof. By the weight reason we have (2) (8) m 

By,®B^~ B^®By,. I 



T\jTj. = Tj.T^ for any r > 1. 



under the isomorphism 



Lemma 3.8. Assume I > k. Identify a two-soliton state with an element of 
Aff{Bi)®Aff{Bk). Then on the set of highest weight elements of Aff{Bi) ® Aff{Bk) , 
the map commutes with the combinatorial R. 



Proof. Con sider Type I case first. We apply T|, on both sides of (3.2). In view of 
Lemma 3.4, we divide the check into 3 cases. 



(i) 5 > 0: 



LHS ^ z"(0 ® z"'-2(d + 1, c, - 1), 
RHS = z"''+"ik) (E) z"'-"-'^{l -b-c+l,c,b-l) 
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(ii) 5 = 0, c> 0: 

LHS = z^'il) ® z'^'-^d + 1, c - 1, 0), 

RHS = z'"'+^(fc) z"'-"-^! - c + 1, c - 1, 0). 

(iii) 5 = c = 0: trivial. 

By Proposition (2), we see the LHS is sent to the RHS by the combinatorial 
R. 

The proof of Type II case gets more involved. Set p = z"^{f) ® z"^ {d,c,b) 
and p' = R{p) = (g) z'"~^(J, 7,/3). Recall integers a, e are defined by 

I = 2e + f , k — 2a + b + c + d. First consider the case when min(a, e) > 0. Then we 
have 



= + 2) ® +l,c,b+ 1). 



(3.4) 

On the other hand, in all cases (i)(ii)(iii) in Proposition 3^ (2) one finds ip < 
k, P + ^ + 6 < I. (Notice that u> < e — a in (iii).) Thus we have 



(3.5) mp') = z"''+"-\ip + 2)®z"'-"-\5+l,j,(3+l). 



One can check ( |3.4| ) is sent to ( p.5D by the combinatorial 7? using Proposition ^ 
(2) again. (Notice that Tt,(p) belongs to the same case as p.) 

Now assume min(a, e) = 0. We divide the cases according to those in Propo- 
sition ^ (2). In all cases we list T[,(p) and ^{p') in the manner ^{p) ^ ^{p'). 
Checking Tif{p) is sent to T'|,(p') by R is left to the reader. 

(i) a > e; We have e = 0, and p' = z"'+^(A:) ® z"^-"{S,c,P) with J = d + Z - 
k — y, l3 — b — y,y — min(^ — k,{b — rf)+), H = k + a+[d — 6)+. 

• bd> 0; The condition implies (35 > Q. 

z'^{l) eg) (d - 1, c, 6 - 1) z'^'+"{k) ® z'^-^-^S - 1, c, /3 - 1). 

• 6 0; We have y ^ 0,6 ^ d + I - k, (3 = 0. 
If c> 0, 

z"(0 ® (d + 1, c - 1, 0) ^ z'"'+^(A;) eg) z'^'^-^S + 1, c - 1, 0). 

If c 0, d < /c - 2, 
2™(0 (g z'"'-i(d + 2, 0, 0) ^ z'^'+^ik) g) z"-^-i((5 + 2, 0, 0). 
If c 0, = fc, 

z'"(0 g) z"'(fc) ^ z"'+^(fc) ® z™--f^(0. 

• & > 0, d = 0, Z - fc < 6; We have y = I - k,6 ^ 0, f3 = b - I + k{> 1). 
Note 6 > 2. 

If c> 0, 

0^(0 ® z™'-\0, c - 1, 6 - 1) ^ z"''+"{k) g) z"'-^~i(0, c - 1, /? - 1). 
If c = 0,/3> 2, 

z™(/) g) z^'-^O, 0, 6 - 2) ^ z"'+^(fc) g) z"-^-i(0, 0, /3 - 2). 
If c==0,/3== 1, 

z™(0 g) z™'-i(0, 0, 6 - 2) z™'+^(fc) g) z'"-^-i(l, 0, 0). 
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• b>0,d = 0,l-k>b;We have y = b,6 = I - k - b, (3 = 0. 

If c> 0, 

® z'"'-i(0, c - 1, 6 - 1) (8) ^'"-f^-i(,5 + 1, c - 1, 0). 
If c = 0,6> 2, 

^'"(0 (g) z"'-i(0, 0, 6 - 2) z'^'+"{k) ® z'^-^-^S + 2, 0, 0). 
If c = 0,6= 1, 

0, 0) z'^'+^ik) ® ^'"-^-1(5 + 2, 0, 0). 

(ii) a < e, > e-a > We have a = 0, and / = 2™'+^(A:)«)z'"--f^(5, c, /3) 
with 5 = d + l — k — e — z,(3=b + e — z,z = mm{b — d + e,l — k — e),H= k. If 
6— d+e < l — k — e, then 2; = 6— d+e, 5> (3 = d, else, ^ = l — k — e, (3 > 6 = d. 

• d > 0; 

_^m-i(j ^ 2) ® z™' (d - 1, c, 6 + 1) z'"'+^(A;) O - 1, c, /? - 1). 

• d = 0,0 0; 

If6-d + e<Z-fc-e, 

2™-i(/ + 2) ® z™'(0, c - 1, 6 + 1) ^ z'^'+^ik) 2™-^-^(5 + 1, c - 1, 0). 

lfb-d + e>l-k-e, 
z'^-\f + 2) ® ^'"'(0, c - 1, 6 + 1) ^ z'^'+^ik) ^'"-^-1(0, c - 1, /3 - 1). 

• d = c = 0; 

If6-d + e</-fc-e, 

z™-i(/ + 2) ® ^'"'(0, 0, fo) z™'+^(fc) (8> z™-f^-i(5 + 2, 0, 0). 

If6-d + e>Z-fc-e,/3>2, 
z™-i(/ + 2) ^"'(0, 0, 6) z'^'+"{k) (8> z™-^-i(0, 0, /3 - 2). 

If6-d + e>Z-fc-e,/3=l, 
^m-i^y + 2) O 2™'(0, 0, 6) z'^'+^ik) ® 2™-^-i(l, 0, 0). 

(iii) a < e,l — k < e — a or e — a < d — b] We have a = 0, and p' = ® 
z'^~"(S. c, /?) with If = k - 2e + 2w,5 = d + I - k - w,P = b + w,w = 
min(/ — k, (2e — d + &)+), Jff = fc + max(e — l + k,d—b — e). Note that w < e, 
and if d = 0, w = Z — A;. 

• d > 0; 

+ 2) ® z""' (d - 1, c, 6 + 1) ^ + 2) ® z'"-^((5 - 1, c, /3 + 1). 

• d = 0, c> 0; 

z'"-i(/ + 2) ® z'"'(0, c - 1, 6 + 1) ^ z'"'+^-i(v5 + 2) ® ^"-^(0, c - 1, /3 + 1). 

• d = c = 0; 

z™-i(/ + 2) (g) z™' (0, 0, 6) ^ 2"'+-^-^ ((^ + 2) ® z™-f^(0, 0, 13). 

I 

Lemma 3.9. Jbr arbitrary l,k G Z>i ('nof necessarily I > k as assumed previ- 
ously), let z^'bi ® z'^ib'i (i = 1,2) be highest weight elements in AfF{Bi) (E)Aff{Bk). 
//T^(z™i6i ® = mz"'-^b2 (8) z'^'^b'^) and wt{bi ® b[) = wt{b2 (S) b'^), then 

z'^'bi (8> z<b[ = 62 O z'^'^b'^. 
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Proof. We first consider Type I case. From Lemma 3.4 we see the inverse image 
of an element of the form z™(Z) z™ {d, c, b) has more than one element only when 
d > 0,b = 0. In this case T^^{z^^{l)^z^"- (d, c, b)) consists of the following elements: 

z"(0®z"'+2(d-l,C,l), 

z"(0®z"'+'(d-l,c+l,0), 
® z™'(d,0,0) ifc = 0. 

It is easy to see the weights are distinct. 

For Type II the cases when r^~i(z'"(/)8)z"''(d,c, b)) has more than one element 
are the following. Using Lemma |3.5| we divide into 6 cases. 

(i) d > 0,6 > 0,a > 0,e = 0: 

z"(0®z'"'+i(d+l,c,6+l). 



z 



m-4-l 



G-2) 



■ z 



(d- l,c,6- 1) 



(ii) d > 0,6 > 0,a 0: 

„ 2) ® z™'+i(d - 1, c, 6 - 1), 
z™+i(/-2)®z™'(d+l,c,6-l). 

(iii) d > 0, 6 = 0, c> 0, a > 0, e = 0: 

®z"'+i(d+ l,c,l), 
^"(0 ®z"'+i(d- l,c+ 1,0). 

(iv) d > 0, 6 = 0^ 



(v) d = 0,6>0,a = 



(vi) d = 0, a > 0, e 



0,e = 






z"(/)® 


z'"'+i(d- 1,1,0), 






z"'+i(d- 2, 0,0) 


if d > 2, 




z'"'+i(0,0,l) ifd 


= 1, 




z"'+i(d+ 1,0,1) 


if a > 0, 




z"'(d, 0,0) ifa = 


0. 


= 0: 








-2)®z"'(0,c+l 


6-1), 


^m+l(j 


-2)®z"'(l,c,6- 


1), 


^m+l(j 


-2)®z'"'(0,0,6) 


if c = 0. 


= 0: 






z"(0® 


z"'+i(l,c,6+l). 




z"(O0 


z"'+i(0,c+l,5 + 


1), 




z"'+i(0,0,6 + 2) 


if c = 0. 



In all cases the weights are distinct. I 

We introduce a degree on the set of highest weight elements AfF{Bi) ® Aif(i?/c). 
First, consider an element p such that Ti^ijj) = p. We define the degree of such p, 
denoted by degp, to be 0. From Lemmas p.4| and we easily see that two-soliton 
states of degree must be of the form z™{l) z™ (fc) for some m and m'. From 
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Propositions |3.1| and |3.2| we also see that any other highest weight elements can 
be transformed to a degree one by applying Ti, for appropriate times. Thus by 
setting degr|,(p) = degp — 1, we can assign a non-negative degree uniquely to all 
the highest weight elements. 



Proof of Theorem 2.7 for Type I and II. Identify a two-soliton state p with 
an element of Aff{Bi) ® AfF{Bk). For r(> k) and sufficiently large t we a re to show 
T^{p) = R{p). Due to the Uq{gn-i) invariance mentioned in Section 3T, it suffices 



to check the assertion when p is a highest weight element. We do this by induction 
on the degree. If degp = 0, the statement can be directly checked by means of 



Lemma |3.3[ In fact the calculation in this case is the same as An'' considered in 
Section |1.2| . Now assume degp > 0. We have 

Here the second equality is assured by the induction hypothesis. Since both T^{p) 



and R{p) have the same weight as p, the proof is finished by Lemma 3.9 



(2) (2) 

3.3. Type III : ^2,^ and D!^_^^ cases. Let us proceed to the case g = 

^2ri I ^n+i- ttlis subsection Bi denotes the J7^(g)-crystal and Bi does the [/^(C^^"*)- 
crystal. 

First we quote some results from |HKOT2|. For g = ^2ri ; ^n+i ^^^'^ ^o^' ^^^y 



I £ Z>i, let uj be the map defined by 
U! : Bi — > B21 

(xi , . . . , Xn ; ■ ■ ■ ; ) 

1-^ (2x1, ■ ■ • , 2a;„, 2x„, . . . , 2xi) for g = 

i-> (2x1, ■ ■ • , 2a;„_i, 2a;„ + Xq, 2x„ + xq, 2x„_i, . . . , 2xi) for g = -D^+i. 
We do not exhibit the ^-dependence in u. The map lu for g,i_i will be denoted as 



Theorem 3.10 ( |HKOT2|l ). Let q ^ A^^,; , dI^I^ . Supposeh(E)c^ c®h under 



the isomorphism Bi ® Bk B^ ® Bi of U^{g)- crystals. Then one has 

(1) uj{b) (Xi oj{c) 1-^ oj{c) ® oj{b) under the isomorphism B21 ^ B2k ^ ^2fe ^ B2 
ofU'q{Cl^^)-crystals. 

(2) H^{b®c)^H^'^'\uj{b)®oj{c)). 

f2l (2) 

For g = ^n+i introduce a map 77 as 



r]-. Bi — > Bi®Bi 

^ 0®0 (ae {!,... ,T}) 

[0] ^ ( only for g = i^i'li). 
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C2) (2) 

Lemma 3.11. For g = and D\_l_^ the following diagram is commutative 



for any L G Z>i . 



B2L ®Bi(g)Bi 



Bi®Bl 



Bi®Bi® B2L 



(2,) (2) ~ ~ 

Proof. For g ^ A\\l (resp. D)^l{) the image r]{Bi) C Bi ® Bi is stable under 
the actions of e1,ff{i ^ 0),eo and /o (resp. ef,ff{i ^ 0, n), eo, e„, /o and /„). 
Moreover these operators generate the identical crystal subgraph within Bi ®Bi as 
they do in uj{Bi) C B2- The assertion follows from this fact and Theorem 3.10| . I 
For Q = A^2n^-^n+i case, denote the map ii : Bi ^ (_Bi)®' specified in Section 
2.3 by ii, and for ci^^ case by ii. Then the following can be checked by a direct 
calculation. 



Lemma 3.12. Letg = A^^j or D^^l^, and\l\ e B^. For any be Bi /orf/^(0„_i) 



l(2 



. p-i^ f«2/('^«-i(&))«)|T| ifl-s{b) is even, 

)) U- j^^^2/(w„-i(6)) ifl-s{h) is odd, 



where s{h) has been defined in Section^ 



Proof of Theorem 2/i_ for Type III. We attribute the proof for q = A'^^I > ^i+i 



to the Cn'' case treated in the previous subsection. Suppose that 



(3.6) 



Z^^bl (g) Z^^b2 ^ Z^^+H2 ® Z^^~% 



under the isomorphis m Aff {Bi)(g)AfF{Bk) ^ Aff{Bk)®)Aff{Bi) of C/q(0„_i)-crystals. 
Then from THEOREM we have 



(3.7) Z«iw„_i(6i) (g) Z^'uJn-l{b2) Z«^+V>-l(fo2) ® Z^'-''LJn-l{bl) 

under the isomorphism AS{B2i) ® Aff{B2k) ^ AS{B2k) ® Aff{B2i) of Uqici^\)- 
crystals. Define e, e', a, a' by 

e — I — s{bi), e' — I — s{bi), a ^ k — 5(62), a' = k — 3(62), 



where s for has been defined in Section 2_^. Put 

Pl = ■ • ■ ® 1 (X) ll{bi)mi ®l® ■ ■ ■ ®l® ife(fo2)m2 ® 1 



e(Si) 



where the notation z;(6i)mi example signifies that vi appearing in ii{bi) = 
Vi (i) ■ ■ ■ ® vi is located on the mi-th component counting from the right end in pi. 
We assume that N — mi, mi — 7712, TO2 ^ I > k. Let t be sufficiently large and 
L > k. Then from Lemmas O, O and |3.11 



U^l «) P2 



we have the commutative diagram: 

Pi <E) uf 

P2 ® ufl 
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Here the upper (lower) ^ denotes the isomorphism of C/^(0)([/q(c4^''))-crystals. 
According to (2.9)-( p!l0 ) the incoming state pi corresponds to 

z^"''-%i z^"'^-%2 e AfF{Bi) (g) AfF{Bk), 

where the RHS is a tensor product of ?7g(0n_i)-crystaL a — {1 — {—l)°')/2 and sim- 
ilarly for e. Thus from (3.6) we are to show that the outgoing state pi corresponds 
to 



(3.8) 



z 



2rrL2—a-i-hi 



Z 



2mi —e — h 



bi e Aff{Bk)®AS(Bi 



To do this we utilize the commutativity of the diagram to explicitly determine pi 
as follows. First we use Lemma 3.12 to find 



P2 = ■■■® '^®^2l{^^7l-l{hl))^,^ ® 1 
(3.9) 

fii ~ 2mi — e, = 2m2 — a. 



1 ® l2fe(Wn-l(&2))^ 



1 



• e Si 



Next we invoke Theorem p. 71 for Cn^ case to derive 



1 



P2 = • • ■ (g) 1 l2/c(w„_l(62))/i2 

(3.10) 

fli ~ fii — 2 min(i, l)t — h, M2 = M2 ~ 2fci + h 



1 «2i(w„_l(6i))^j (g) 1 



(«.) 



®2N 



where t he ph ase shift is h in ( |3.6| ) because of Theorem 3.10 (2). Finally by using 
Lemma 3.12 again we obtain 



(3.11) 



Pl = • • • (g) 1 l;(62)m2 g) 1 (g • • • (g 1 ® i/c(&l)mi ® 1 ' 



Ml + e' 
™i = , "i2 



Ai2 + a 



According to (2.9) and (2.10) this corresponds to 



^2(m2 + kt)-a'^^ ^ ^2(mi+min(L,/)t)-e'^^ 



Combining (|3.9|),(3.1C) and (3.11), we find that this is indeed equal to (3 



4. Examples 

Let us present examples of soliton scattering. Given a state p of an automaton, 
to find the time evolution Tr{p) involves a series of calculations oiBr®Bi ~ Bi®Br- 
For r finite this is a tedious task in general but can be done according to the insertion 
scheme in [ HKOTl , HKOT2| . The presented examples are actually so generated 
by a computer by setting r = 12. However if r is large enough so that T = Tr, the 
result in [ HKT2| ] enables one to compute T(p) easily in an analogous manner to 
the An'' case where the simple algorithm (box-ball rule) is known | |TS| , |TNS |. 
In all the examples below T = T12 is valid. 

First we give examples of two-soliton scattering described by AfF{B^)(gAfF{B3) — 
Aff{B3)^Aff{B5). 

Example 4.1. g = a''^\ 

: • • • 11^ 44322 ^11111^433^111111111111111111111111111111 • • • 

1 : • • • 1111111^ 44322 ^11^3^111111111111111111111111111 • • • 
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111111111111^ 44322 ^1^3^111111111111111111111111 
11111111111111111^ 4432^433 ^111111111111111111111 
llllllllllllllllllllli p2 ^ 44^4 pillllllllllllllll 
111111111111111111111111132 ^14^4 ^433111111111111 
1111111111111111111111111111322^ YY44 ^4331111111 



Here the markers specify the positions of initial soHtons under the interaction- 
free propagation. One can read off the phase shift from the deviation of the outgoing 
sohtons from the corresponding markers. In terms of the soHton labels, the above 
scattering is described by the combinatorial R of [/g(yl2^-')-crystal: 

£"(2, 1, 2) z"'(0, 2, 1) h-> z'"'+3(2, 1, 0) z""%0, 2, 3). 



Example 4.2. g = A)' 



(2) 



11^4222|111111|222|11111111111111111111111111111111111111 
1111111 ^4222 ^111^22^11111111111111111111111111111111111 



111111111111 ^4222 ^11^2^11111111111111111111111111111111 



11111111111111111^4222^11111111111111111111111111111 
llllllllllllllllllllli p42 ^ lip ^llllllllllllllllllllllllll 
lllllllllllllllllllllllllll ^lil ^ ^ llllllllllllllllllllll 
lllllllllllllllllllllllllllllll^p24|lllllllllllllllll 
11111111111111111111111111111111112222222^111111111111 



1111111111111111111111111111111111111^22^ ^2224^ 1111111 
1111111111111111111111111111111111111111^22^111 ^224 ^11 

Compare this with the combinatorial R of Uq{A\^ ^)-crystal : 

^"(3, 0, 0, 1, 0, 1) (g> z™' (0, 0, 0, 0, 0, 3) 

t->z"'+°(3, 0, 0, 0, 0, 0) (E> z"'-"(0, 0, 0, 1, 0, 4). 

The f ollow ing example starts with the formally identical initial condition as 
Example L2 but under a different g. The outgoing solitons are diffe rent. Moreover, 
separation into two final solitons is much quicker in Example 4.3. It has a larger 
phase shift, which corresponds to the greater value of the energy functions of the 
combinatorial R. 



Example 4.3. 



— ^8 • 



11^4222|111111|222|1111111111111111111111111111111 
1111111 ^4222 ^111^2^1111111111111111111111111111 
111111111111 ^4222 ^11^2^1111111111111111111111111 
lllllllllllllllllp4222^1111111111111111111111 
llllllllllllllllllllllp42|lglllllllllllllllllll 
111111111111111111111111124 p ^ lll ^ i pllllllllllllll 
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6 : • • • 1111111111111111111111111111242^ YYll plllllllll • • • 

7 : • • • 1111111111111111111111111111111242^ 11 Yllll pilll • ■ • 

(2'] 

Compare this with the combinatorial R of C/g(Ag )-crystal : 
z™(3, 0, 0, 1, 0, 1) £"'(0, 0, 0, 0, 0, 3) 
^z'^'+^{\, 0, 0, 1, 0, 1) ® z"'^'^(0, 0, 0, 0, 0, 1). 

Not e that A^^^ is Type III, so the counting of the phase is doubled as in ( ^.9|) 
- ( 2.1C ). In particular when G i?i locates in front of solitons, there is a special 
feature ( 2.10 ). The following is such an example. 



Example 4.4. q — A 



(2) 



11 ^4333 ^11111^042^11111111111111111111111 
1111111 ^4333 ^11^042^11111111111111111111 
llllllllllllp4333|l|042|lllllllllllllllll 
lllllllllllllllllp431|^2111111111111 
1111111111111111111124^ 31.110,4 3321111111 



111111111111111111111112431^ 11^0 ^433211- 

Compare this with the combinatorial R of [/g(j4g^-')-crystal : 
2™(0, 3, 1, 0, 0, 1) £"'(1, 0, 0, 1, 0, 0) 
t^z'"'+9(0, 1, 1, 0, 0, 1) z"'-^(l, 2, 0, 1, 0, 0). 



Example 4.5. q = b''1\ 



1111111^ 44433 ^11111^440^1111111111111111111111111111111 
111111111111^ 44433 ^111^440^1111111111111111111111111111 
11111111111111111^ 44433 ^1^440^1111111111111111111111111 
1111111111111111111111^443^40|1111111111111111111111 
111111111111111111111111111^ 444^4p p311111111111111111 
11111111111111111111111111111144^^1|4|033111111111111 
111111111111111111111111111111111444|1^1|144|0331111111 
111111111111111111111111111111111111444^ 1iy 114^ 03311 

Compare this with the combinatorial R of C/g(i?3^'')-crystal : 
z™(0, 2, 3, 0, 0, 0, 0) ® z"' (0, 0, 0, 1, 2, 0, 0) 
H^z"'+^(0, 0, 3, 0, 0, 0, 0) ® z""3(0, 2, 0, 1, 2, 0, 0). 



Example 4.6. q = C, 



(1) 



• • 11^ 12341 ^111111.243.111111111111111111111111111 • • ■ 
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1 : • • • llllll Y234Y llip43^111111111111111111111111 • • • 

2 : • • • 11111111111 Y234Y ip43^111111111111111111111 • • • 
3 : • ■ • 11111111111111111^ 12311^44 ^11111111111111111 • • • 
4 : • • • 1111111111111111111111|311^443111111111111 ■ • • 
5 : • ■ • 111111111111111111111111131^ Yll^ p4431111111 • ■ • 
6 : • • • llllllllllllllllllllllllllli311^ YY12 p44311 • • • 

Compare this with the combinatorial R of J7g(C3^'')-crystal : 

2™(0, 0, 0, 1, 1, 1) z""' (0, 1, 0, 1, 0, 1) 
^^m'+4(Q^ 0, 0, 0, 1, 0) O ^;™-^(0, 1, 0, 2, 0, 2). 

Example 4.7. g = D^^K 

: ■ ■ ■ 11^ 35322 ^111111^5^111111111111111111111111111111111111111 
1 : • • • 1111111 ^5322 ^11ip25^111111111111111111111111111111111111 

2 : • • • 111111111111 ^322 ^11^5^111111111111111111111111111111111 

3 : • • • 11111111111111111 ^322^25 ^111111111111111111111111111111 

4 : • • • 1111111111111111111111 ^3 ^ 11 ^ 5 ^111111111111111111111111111 

5 : • • • 111111111111111111111111111^ 4^114 ^ 5 ^51111111111111111111111 

6 : • • • 1111111111111111111111111111114 ^2 ^ 245 ^1111111111111111 

7 : ■ ■ ■ 1111111111111111111111111111111114 ^2^2245 ^5111111111111 

8 : • • • 1111111111111111111111111111111111114^21^ Y2245 ^51111111 

9 : • • • 11111111111111111111111111111111111111142211111.12245.511 



Compare this with the combinatorial R of Uq{D^ ■^)-crystal : 

^™(2, 1, 0, 0, 1, 0, 1, 0) (g) z""' (0, 0, 0, 0, 1, 0, 0, 2) 
h^^™'+^(2, 0, 1, 0, 0, 0, 0, 0) O ^™-^(0, 0, 0, 0, 2, 1, 0, 2). 

Example 4.8. g = D^^\ 

: • • • 11 ^0322 ^111111^101^111111111111111111111111111 • • • 

1 : • • • 1111111 ^0322 ^111^101^111111111111111111111111 • • • 

2 : • • • llllllllllllp0322|ll^lllllllllllllllllllll • • • 

3 : • • • 11111111111111111 ^0321 ^ 102 ^111111111111111111 • • • 

4 : • • • 111111111111111111111^031|l|0p21111111111111 • • • 

5 : • • • 111111111111111111111111203^ 1 ^ 11 Y p2211111111 • • • 

6 : • • • 1111111111111111111111111112031.1.11.111.022111 ■ • • 
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Compare this with the combinatorial R of Uq{D^^^ )-CT-ystal : 
z"(2, 1, 0, 1, 0, 0, 1) (g> z™' (0, 0, 0, 1, 0, 0, 0) 
h^z'"'+®(0, 1, 0, 1, 0, 0, 1) (E> z"-«(2, 0, 0, 1, 0, 0, 0). 

(2) 

Note that I?^ is Type III, so the counting of the phase is doubled as in 
Examples O and 



Finally we present two comparable figures showing the soliton nature, namely, 
independence of the order of the collisions. 



Example 4.9. ^ = -84 



(1) 



llp22|lllgllllllllQlllllllllllllllllllllllllllll 
lllllp22|ll^lllllllQllllllllllllllllllllllllllll 
llllllllp2|l^llllll^llllllllllllllllllllllllll 
lllllllllllp3l|3^1111|0|llllllllllllllllllllllllll 
1111111111111^ 31^ 3311^111111111111111111111111 
11111111111111123^ 1^ 31^011111111111111111111111 
11111111111111111231^^1^3011111111111111111111 
1111111111111111111231 ^1^ 1133011111111111111111 
111111111111111111111231,31,1.11,1133011111111111111 
11111111111111111111111233^1^111^133011111111111 
1111111111111111111111111323^11^111^133011111111 
1111111111111111111111111131^3^11^111^111^133011111 
111111111111111111111111111311^^1111^111^133011 

llp2|llllllllllglllQlllllllllllllllllllllllllll 
lllllp22|lllllllllgl 1^11111111111111111111111111 
llllllllp2|llllllllgl^llllllllllllllllllllllll 
Illllllllllp22|lllllll|31j3|011111111111111111111111 
11111111111111^22^111113^0111111111111111111111 
llllllllllllllllip22^11131^01111111111111111111 
11111111111111111111^22^131^1^11^011111111111111111 
11111111111111111111111^33^3^1^0111111111111111 
11111111111111111111111113^ 1^ 31^01111111111111 
11111111111111111111111111311^3|3^011111111111 
111111111111111111111111111311^^123^3011111111 
1111111111111111111111111111311^11 ^3^ 33011111 
11111111111111111111111111111311^1112^ 3^ 1133011 



The both scatterings are described by 

z'"(l, 0, 0, 0, 0, 0, 2) ® z"'(0, 0, 0, 0, 0, 2, 0) ® z""(0, 0, 0, 1, 0, 0, 0) 



(0, 1, 0, 0, 0, 0, 0) z™ +^(0, 0, 0, 0, 0, 1, 1) ® z™'^(0, 0, 0, 1, 0, 2, 0) 
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under the combinatorial R of C/g(S^^')-crystal : Aff{B3) ® Aff{B2) AfF{Bi) ^ 
Ati{Bi) (g) Aff(B2) Aff(B3). 

5. Discussion 

So far we have only dealt with the homogeneous tensor product • ■ • ® -Bi ® -Bi 
• • • with a ferromagnetic boundary condition as the automaton states. It is natural 
to consider a generalization such that the states ■ ■ - ^hj® hj+i are taken from 
the inhomogeneous tensor product ■ • ■<E)Bmj ^B^j+i®- ■ ■ , where m^-'s are arbitrary 
positive integers. See [HHIKTT| for A\^^ case. Imposing the boundary condition 



= Urrij (see (2.1)) for all the distant j, the time evolutions can be defined in an 



analogous manner to Section 2.1. Let us discuss solitons and their scattering rule 
in the inhomogeneous case. 

First, in order to define solitons, we actually consider the tensor product of the 
form: 

■ ■ ■ (g) Bi (g) Bi (g) ® • • • B^„ (g) Bi (g) Bi (g ■ ■ ■ . 

Here the left region ■ ■ -(gBi^Bi and the right region Bi(gBi®- ■ ■ are supplemented 
as the detectors of incoming and outgoing solitons. In those regions one can assign 
elements in Aff(i3/j^) (g) • ■ • (gAff(Bi.) to any asymptotic (enough separated) i-soliton 



states in the same way as ( 2.11 ) . One caution necessary here is how to specify the 



phase of solitons in the left region. To count the "distance" "f{t) from the right 



(see after (2/7)), we understand that each component Bm^ in the middle region 
is effectively occupying mj slots. In this way, to a soliton of length I in either 
asymptotic region, one can assign the data z'^b with b from the C/^(g„_i)-crystal 
Bi. 

Now we are ready to study the collisions of solitons which take place in the 
middle inhomogeneous region in general. Consider the incoming two-soliton state 
represented by 

z'-'^bi ig z^^b2 e AS{Bi) (g AfF{Bk), (l > k) 

in the above sense, and let z'^^+%®z'^^-% e AS{Bk)®AS{Bi) with 5 = H{bi®b2) 
be the image under the combinatorial R : AS{Bi) ® AS{Bk) ^ AS{Bk) ® Aff{Bi). 
Then we claim that after sufficiently many time steps under the evolution Tr{r > fc), 
the state again becomes the two-soliton state corresponding to 

^72+<5-A.^^ (g z-'^-^-^'bi E AS{Bk) eg Aff{Bi), 

N 



Comparing this with Theorem 2.7, we see that the effect of inhomogeneity only 



enters the additional phase shifts A/ and A^. In fact the claim reduces to Theorem 



2.7 when Vm^ = 1. To justify the claim let us first note that for any b E Bi and 



L, L' G Z>o, the relation 

HL+Xi) I ,(g)L 



(g) li(6) (8)|JJ g) Umi g) • ■ • ® Mm„ 



)|T| (gH{b)(g[T] 

sidering the evolution of the soliton ii{b) under Tmi ■ ■ ■ and applying Lemma 



is valid, where Xi = iiiiii(z, mj)(= X^fci ™j ~ A^). This can be shown by con- 
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p.5[^(2). It means that the sohton ii{b) gains the phase shift L' - 

= when penetrating the inhomogeneous region 
We depict the resuft as 



X,, 



-{L' + 



(5.1) liiV) MjTii ' ' ' Wmjv ^ Umi'''UmN 

for some t 1, where ... denotes the tensor products of [T] G Bi. 

Now we proceed to the consideration on the incoming state in the claim. In 
the above notation it looks as 

Pi = n{bl) Ik{b2) Umi •••Wm„ i 

where we have suppressed the initial phase 71 , 72 ■ We shall keep track of their shifts 
only. The collision may take place in the middle region. However by means of ( ^.2|) 
one can push it away to the right as far as one wishes. Namely pi ~ p'l, where 

p'l = nibi) 1^(62) ■■■Um^ 

Pushing it away enough, one can let two solitons fin ish the collision within the 
enlarged left region. To such a collision Theorem 2/? is applied, leading to 

s -s 

for some ti ^ 1, where the phase shift is given hy 6 = H{hi ® 62). Now we can 
utilize to let the two solitons penetrate the inhomogeneous region. This can 
be done individually for each soliton because the two solitons in the above may be 
assumed to be enough separated and the right one is faster than 1^(62) under 

Tr with r > k. Thus for some ^2 ^ we have 

Tr^ip'l) = ""mi •••Mmjv «A:(&2) M^l)---- , 

where the solitons have gained the extra phase shifts —A; and — Afc. Fi nally we 
can pull the region Um^ 
result we have T^^{p[) z 



back to the original position again by ( |2.2| ). As the 



P2, where 



-S-A, 



P2 = Wmi • ■ -Umjv Ikih) , 

and this tensor product has the same ordering of the components Bi 's and B„i- 's 
as pi. Since Tr commutes with the isomorphism ~, we conclude Tjf^ (pi) = p2. This 
completes the proof of the claim. 

Example 5.1. Consider the q = Cg^'' automaton with the inhomogeneous re- 



gion TOl 
length / : 



= 3, m2 = 
3 and k 



2 {N = 2). The following is the collision of two solitons with 
= 1 in the inhomogeneous region under the time evolution T12: 
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1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . Ill . 1^. 2 . 1 ■ ^ 1 ■ 1 ■ 3 ^ ■ 3 ■2-1-1-1-1-1 
1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . Ill . 11 2 . 1 . 1 • 1 ■ ^ 1 ■ 1 ■ 3 ^ ■ 3 • 2 • 1 • 1 

Here^_stands for ®, and 111 := (3,0,0,0,0,0) £ Bg, 112 (2,1,0,0,0,0) £ 
^3, 132 (1,0,0,1,1,0) e ^3,12 (1,0,0,0,1,0) G B2, etc. On the other hand, 
one has 

z™(0, 1, 1, 1) 8) (1, 0, 0, 0) ^ z™'+2(0, 0, 0, 1) z"'-2(l, 1, 1, 0) 

under the combinatorial R of ;7,(C^^^)-crystal : AfF{B-i) (g) Aff{Bi) ^ AfF{Bi) 
Aff(_B3). Since A3 = and Ai = 3, the phase shifts of the larger and smaller 
solitons should be —2 — A3 = —2 and 2 — Ai = —1, respectively. These values 
indeed agree with the deviation of the outgoing solitons from the corresponding 
markers. 
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Appendix A. Crystals and the combinatorial R 

We review fundamental facts on crystals. It is assumed that our crystals are 
crystal bases of irreducible finite-dimensional representations of a quantum affine 



algebra U'{q). Our basic reference is [KMN 



Let P be the weight lattice, {ai}o<i<n the simple roots, and {Ai}o<i<„ the 
fundamental weights of g. Let denote the classical part of A^. See Q for 
details. A crystal B is a finite set with weight decomposition B = U\^pB\. The 
Kashiwara operators Ci, fi (i = 0, 1, • • • , n) act on B as 

h-.Bx^ Bx+a. U {0}, /. : Sa ^ Ba-o. U {0}. 

In particular, these operators are nilpotent. By definition, we have fib = b' if and 
only if 6 = iib'. Drawing 6 &' in such a case, B is endowed with the structure of 
colored oriented graph called crystal graph. For any b G B, set 

e,(6) = max{TO > | i'^b 7^ 0}, ^^{b) = max{m > | /ffc ^ 0}. 

Then we have the weight wt6 of 6 by wtb = X]r=o('^«(^) ~ £i(^))^i- 

For two crystals B and B', one can define the tensor product B®B' — {b®b' \ 
b e B,b' e B'}. The operators ii, fi act on B ® B' by 



ei{b®b') 
h{b®b') 



eib®b' if Lpi{b) > Si{b') 
b®e.ib' if (pi{b) < e.i{b'), 

~f,b®b' \i ^,[b) > e,{b') 
b®j,b' if ^,(6) <£,(&') • 

Here ® &' and b®Q should be understood as 0. For crystals we are considering, 
there exists a unique isomorphism B ® B' ^ B' (i) B, i.e. a unique map which 
commutes with the action of Kashiwara operators. In particular, it preserves the 
weight. 
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For a crystal B we define its affinization Aff{B) = {z'^b | c? £ Z, 6 G B}. 
z is called the spectral parameter. AfF{B) also admits the crystal structure by 
e,-z% = z'^+^^°{e,b)Ji-z% = z'^'^^" {f,b) . The crystal isomorphism ^ B'(g)B 

is hfted up to a map AfF{B) » AfF{B') ^ AfF{B') (g) AfF{B) called the combinatorial 
R-matrix (combinatorial R for short). It has the following form: 

R : Aff{B) (g) AfF{B') — > AS(B') ® AS{B) 

z% ® z^'b' ^ zd'+H{b®b')y ^ ^d-H(b»b')i^ 

where b^b'i-^b'^b under the isomorphism B ® B' ^ B' ® B . H{b ® b') is called 
the energy function and determined up to a global additive constant by 

r H{b ®b') + l if ^ = 0, (^o(fo) > eo{b'), Mb') > eo(6), 
H{e,{b ® 6')) = < 6') " 1 if * = 0, Mb) < £oib'), Mb') < £o(b), 

[ H{b (g) b') otherwise . 

Now we can state 

Proposition A.l (Yang-Baxter equation). The following equation holds on 
Aff{Bi) ^ Aff{Bi,) ® AS{Bin): 

{R R){R 1) = (1 i?)(i? ® 1)(1 R). 

Suppose I > k. In the main text we fix the values of the energy function H on 

Bi ® Bk so that 

(A.l) Hiui Ufc) — 2<;fc for any g, 



regardless of the rank of g. See (p.5|) for ^. The values given in PROPOSITIONS 3.1 



and |3.2| are so normalized. In fact max{iJ(6 ® &') \ b € Bi,b' ^ B^} = 2<;k holds, 
and as for the minimum we have 

min{H{b®b') \ b e Bi,b' E Bk} = fc, 0, 2fc, 0, fc, 0, 2A; 

Appendix B. Crystal Bi 

Let us present the parameterization of the crystal Bi as a set for each q. 
3^A^^: 



ri+l 



Bi = {(xi, . . . , x^M e Z"+i|a;, > 0, ^ x, = I}. 

i=l 

For i?i we use a simpler notation 

(xi, . . . , a;„+i) = [T] a Xi = 1, others = 0. 

— ^2n-l • 

n 

Bi = {(xi, . . . ,a;„,x„, . . . ,xi) e Z^"|xi,Xi > 0, ^(a:^i + Xj) = 0- 

1=1 

For i?i we use a simpler notation 

I [T] if Xi = 1, others = 0, 



(B.l) (xi,...,x„,x„,...,xi) 



[ [T] if Xi = 1, others = 0. 
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Bi = {(Xi, . . . ,Xn,Xnj ■ ■ ■ ,Xi) G Z^" | Xi, Xi > 0, '^{Xi + Xi) < 1} 

i=l 

For Bi we use a simpler notation 

{ [T] if Xi = 1, others = 0, 
[T] if Xi = 1, others = 0, 
if Xi = 0, Xi = for all i. 

n - r(i) ■ 



B/ = |(xi, . . . ,x„,xo,x„, . . . ,xi) e Z^" X {0,1} 
For Bi we use a simpler notation 



xo = or 1, Xi, Xi > 0, 

2^0 + J27=li^i +Xi) =1 



(Xi , . . . , X^ , Xq , Xj^ , . . . , Xi ) 



— ^ J if ^4 = Ij others = 0, 
[ [T] if Xi = 1, others = 0. 



= : 



= < (xi, . . . ,x„,x„, . . . ,xi) e Z- 



2n 



For Bi we use the notation ( |B.l| ) 
= Z?.i^': 



5/ = <^ (xi, . . . ,x„,x„, . . . ,xi) e Z- 



2n 



x„ = or x„ = 0, Xi, Xj > 0, 

J2"^l{Xi+X^) = I 



For Bi we use the notation ( |B.]| ) 



5; = |(xi, . . . ,X„,Xo,X„, . . . ,xi) e Z^" X {0,1} 
For Bi we use a simpler notation 



Xq = or 1, Xi, Xi > 0, 

+ z]r=i(^* + Xj) < z 



Bi 9 (xi,...,x„,xo,x„,...,xi) 



[T] if Xi = 1, others = 0, 
[T] if Xi = 1, others ~ 0, 

if Xi = 0, Xi = for all i. 
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Appendix C. Crystal structure of 

We review the crystal structure of _Bt] for Type I and II, i.e., for q — A^^_i, Bn^ , D 
Define an index set J and a linear order ^ on J by 

J={1<2< <n<Q-^n< <2^T} for0 = B,W, 

= {l-<2-<----<n-<7l^---^2^T} for0 = 4^^^^, , C^^) . 

For g of Type I or II, let g be the classical part specified in the beginning of Section 
1.3. As a set is given by 

- B{K2) © B{0) for Type I (g = i , , 
= B(A2) for Type II (g = &^^), 
5(0) = {<^}, 



B{K2 



'a 

(1) a-<(3 if g = C„, 

a^/3or (a,/3) = (0,0) ifg_=B„, 
a ^ /3 or (a, /3) = (n, n) if g = Dn 

(2) (a,/?)^(l,T). 



a, P € J satisfying conditions 
(1) and (2) below 



The labeling of the elements of B{A2) is taken from |KN]. Needless to say, the 
actions of the crystal operators e^, /i(i = 1, 2, • • • , n) on B{A2) agree with those in 



KN[ |. On B{0) they are defined by 

Szcl) ^ fi<j> = (i = 1,2, • • • ,n). 

Therefore we are left to show the actions of eo, fa, which are given as follows: 
For Type I, 

a ^ c: 

'2' 



a 



1/ \a 

a 
P. 



if a 7^ 1,2, 
if a 7^ 2, 2, 
in the other cases. 



For Type II, 



if 1, 
in the other cases. 



As usual, b — ^ b' means fib = b' <==^ b = Cib'. 
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